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Abstract 


developments in high energy lasers, adaptive 
\ 0 P^cs, and atmospheric transmission bring laser oro- 
V^^on moch closer to realization Perhaps more* 
impiortant, the need to transport much greater tonnaaes 

^d fSr t mX? 1 -t^ 0mmerClal purposes ’ s Pace Station Freedom 
"^v tary P ur POses is now clear. A part (e a 
JjjJf H tbe sp M® station supplies) of this traffic 
could be orbited in small packages. Accordingly a 
workshop on thfs possibility was held at Livermore 

1986 and th,s wer uans 

/ J t( ’ 4 ( ~ *'-**£ 4 - ^ " 

^ * 8 P a P® F . -proposes a reference vehicle for studv 

Tee? TlluL? ? ay,0ad ? nd so,id Prop.ir.iJ e.g. 

icej. a suitable laser pulse is proposed for usina a 
efficiently*^^ Detonation wave to Produce thrust 9 

llk ? ly t t iat a min i"ium system (10 Mw C02 laser 
tiin M d Tk^ 1rr ° r - could be constructed for about 

J hl ? minimum system could launch payloads of 
about 13 kg to a 400 km orbit every 10 minutes, 
he annual launch capability would be about 683 tonnes 

order ^f 6 duty factor - Laser propul sion would be an 
the d„? f J ag Jl nude ch eaper than chemical rockets if 
the duty factor was 20% (10,000 launches/yr. ) and 
launches beyond that would be even cheaper. 

The chief problem which needs to be addressed before 
these possibilities could be realized is the desiqn of 
a propellant to turn laser energy into thrust effi- 
ciently and to withstand the launch environment. f 


INTRODUCTION 


The key cost which determines the magnitude of realistic doss1M11m»c 
in space is th e cost of transportation to low Earth orbiMLEO) One of 
the great disappointments In the utilization of spSce is that in the 29 


years since Sputnik this key cost has not declined. 

One opportunity for dramatic Ijprov^jt can 

energy from a laser on th ® Jj r °V. . ible t Q transfer energies to 

easily vaporize any material and it is possib e to trans ^ ev ^ porate d 

the vapor which are ! a [9 e ° the vehicle, and the kinetic energy 

of launching a payload every few minutes. 

uhon * Mc - v «. tpm was first proposed fifteen years ago (Refs. 1 and 2) four 
“ojSr extrapolations of existing technology were required for its 
implementation: 

1. Laser average powers had to be extended by several orders of magnitude. 

2. Atmospheric transmission problems needed to be explored. 

3. Collimating mirrors larger than conventional dimensions needed to be 
developed. 

loSo) Md wi?S thrust^vectors consldorably off the laser beam axis 
needed to be developed. 

'v-' n* 1 « ior 

program was undertaken in the early seventies. 

SL’xis £ c t«r» c ptn n9 

megawatts through the atmosphere are being addressea„ n „ e other 

SDi^^r^rtS^^IreV^su s^ Nations tha^these^eam^can be^ 

MS! Ms'ISdSTSSibli to build very large mirrors, e.g., 
the 10 meter Keck telescope. 

The development of thruster technology has not been * P“«"Jf 
changing MMaSIm! jSiS’ST'Jh Smaller laser power). 


♦Strategic Defense Initiative 
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Phases 1 & 2 


Phase 3 


Pay toad 





Laser beam 



s:::;r. tss: * - ’"' , * -■ - - •*• ^assssa; 

face as rapidly as possible. The plasma formed will then cKiaIh tha 
pellant from the LSD phase which fallows. Shleld the pr °' 


I would like to propose the modification of this two pulse system illus- 
trated in Fig. 1A. Here the pulse is divided into three phases with the 
addition of an "ignition spike." The evaporation energy is typically 
one to two orders of magnitude smaller than the energy needed to dnve 
the LSD The soecific energy deposited in an element of gas by the LSD 
U proSirt!Sn.? P " the 2/3 power of the ratio of the flux to the ens ty 
(Ref 3) Thus the deposited energy can be controlled either by control- 
llngthe density or the flux during the LSD phase. It is clear y inore 
economical to control the density distribution which depends only on the 
flux as a function of time during the evaporation phase. It is therefore 
proposed that the flux during the evaporation be an adjustable function 
of time chosen to produce a density distribution designed to achieve the 
desired specific impulse with maximum efficiency. 

The ignition phase consists of a spike, perhaps a gain switched spike, 
intended to ionize the vapor close to the propellant surface sufficiently 
to make it opaque to the laser radiation. It is important that this pro- 
cess be accomplished rapidly to avoid evaporating too much propel ant at 
this staqe. It will be seen that the highest attainable flux will lead 
to minimum evaporation (see reference 4). It is therefore proposed 
that the flux in the ignition spike be as large as can be delivered. The 
limitation will probably be imposed by surface breakdown in the laser. 
Required duration of the ignition spike will depend on the nature of the 
propellant e.g. for lithium Hyde calculated that, at a flux of 25 Mw/cm^2, 
thetime required was 7.5 nsec and thus a fluence of less than .2 joules/ 
cm A 2 to achieve unit optical depth. 

The third LSD phase involves most of the energy and the cost of creating 
it controls the system cost. The objective of this pulse shape design has 
been to prepare the vapor to efficiently utilize whatever pulse shape 
minimizes this laser cost. 

ESTIMATES OF THE CAPABILITY OF THE REFERENCE SYSTEM 

It will be useful to start by attempting an estimate of the losses which 
are foreseen for this process. 

First there are the losses in the laser itself. Assuming electrically 

driven lasers, the "wall plug" efficiency EL gives * h « 

in the collimated laser beam to the power drawn from th J 

For example for a 10 micron free electron laser Briggs ( LLNL) informed us 

that an appropriate EL would be about 20%. For the C02 laser Daughtery 

(AVCO) suggested 16%. 

Second there are losses in the atmosphere. Among these are scattering 
due to Thermal Blooming, Stimulated Raman Scattering, and Atmospheric 
Turbulence, and there is absorption chiefly due to water vapor. e 
assumed that at the intensities to be used (for the sample trajectory 
10 A 5 w/sq cm at the collimator) and noting that we have a cooperative 
taraet/that the beam would be essentially diffraction limited. Starting 
fr^nZtein top about 10,000 ft high we took the trans.,ss,on t rough 
the remainder of the atmosphere, EA, to be .9/cos(TH) where TH is the 
zenith angle. 
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The third set of losses occurs in converting the laser energy arriving at 
the vehicle into thrust. An ideal thruster would convert the laser * 
energy arriving at the vehicle into kinetic energy MD0T*VJ A 2/2 where MDOT 

J2«2 e *!! r ! p ? 11 !E t !J se 2 1 per , second and VJ 1s the de s* 9 ned jet velocity. 

(Note that in the double pulse system, the propellant mass and the pro- 
pulsive energy can be chosen independently; VJ can be chosen to optimize 

eff[ciln?r f ETH an include° SSeS ^ COnversion whtch result in * thruster 

A. The latent heat of the evaporated propellant. 

B. Chemical or internal molecular energy remaining In the jet following 
the one d i me nsi onal expansion. These losses will be reduced by the use 

l0nges - du r a J i0 ? P ulses wh1ch still allow essentially one dimen- 
sional expansion ( 1 microsecond for a 1 meter dla vehicle). 

C. Losses due to non-homogeneities in gas velocity in the jet. 

These losses are estimated in reference 5 by Rod Hyde for 

lithium as a propellant. (Lithium was chosen simply for ease of calcula- 
ti on • j 

in t 5 e , value of ETH are presently the leading uncertain- 
th*. wn!tch efflC1enCy ° f aser P r °P uls1on - At the present state of the art 

(which was «d\ aW in°Re7 aS l 0 ) n t0 Ch< "' 9e 9 “' SS that ETH ' rou,,, be about 401 

finally the kinetic and potential energy in the payload is of course 
smaller than the kinetic energy in the propulsive jet integrated over the 
trajectory. We will call thematic, of these energies tSe trajectory 
efficiency, ETR. For the sample trajectory ETR was 27%. y 

ESTIMATE OF THE MASS LAUNCHING CAPABILITY, Mo, OF A PULSED LASER 

The ra nge to orbit, D, will be dominated by the acceleration during the 

vnnT 1ty V p ° rt l ons ? f the trajectory. If we take this accelera- 
tion, VDOT, as constant, using the final acceleration and taking VJ = V 


VD0T*Mo = MD0T*VJ = 2*ETH*P‘/V, 
we get 

n = n Mo*V*3 
u 2*VD0T wra*p’ 


(1) 


( 2 ) 


Mn . . . where P’ is the average laser power at the vehicle, and 

Mo is the mass launched. 

The radius, rv of a vehicle base which can be Illuminated with a flux + 
with peak power PP from the laser Is 9 

rv * (PP/pl* 4 )*.5. 

“? top ” d1str,butl0n of intensity. Actually 0 will 

be limited by diffraction to approximately 
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Dmax = rm*rv/.3*x 


(3) 


where rm is the radius of the collimating mirror, A is the 
laser wavelength and the constant .3 is chosen to correspond to a value 
of rv ha If way 9 to the first dark ring. For A = 10 microns, rv = 50 cm, 
and rm = 5 m we get Dmax = 833 km. 


Setting D = Dmax (eq. 2&3) we get 


Mo = 


7. 5*rm*ETH*P , *(PP/ ♦ T-5 

X * V*3 


(4) 


Allowinq for atmospheric absorption as discussed above we take for our 
model 10 MW laser, P’ = 8 MW, PP = P’*10 A 4, ♦ * 1CT7 watts/cm 2 and the 
above values of rm, rv, ETH, and 1, the mass which can be accelerated to 
orbital velocity is 18 kg. The agreement of this estimate with the 
sample trajectory result calculated below (13.79 kg.) is as good as 
might be expected since in that case the acceleration was not constant. 

A SAMPLE TRAJECTORY* 

When we consider the practical applications of laser propulsion an im- 
portant consideration is the minimum** scale of an initial trial. 

To illuminate the choice for this minimum scale we will attempt to cal- 
culate the payload which can be launched with a 10 MW laser, making the 
guess that this will be within one order of magnitude of the practical 
minimum. It was assumed that the pulse duration was 10^-4) J;he time 
between pulses so that the flux on the 1 sq meter vehicle base would be 
10 MW/sq cm. Note that the mass which can be launched varies inversely 
with the square root of the minimum flux sufficient to sustain an effi- 
cient LSD. The achievement LSDs at low flux will be one of the most 
important objectives of propellant research. 

After several trials it was found possible to launch 13.79 kg into a 411 
km circular orbit making the assumptions listed in Table 1. 

* In Ref. 5 Jordin Kare gives a more complete modelling 
of the laser launching. Close agreement between his results and those 
presented here provides some confidence that the remaining bugs are not 
too important. 

** We will not consider here the utility of small satellites other than 
to note that Freeman Dyson* proposes that satellites as small as 1 kg 
would be useful for space science purposes. 


I 




♦Dyson F., see his March 26, 1986 talk at Analog Devices, 
Norwood, MA. 
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TABLE 1 


I. 


ASSUMPTIONS FOR SAMPLE 10 MW LASER LAUNCH MODELLING 

A. Initial mass (propellant + payload) * 120 kg. 

B. The base area of the propellant = l S q. meter 
c. The coefficient of atmospheric drag = .4 

sphere!*) S aSSUmeS ^ th * Vehicle Wl11 be streamll ‘ned as well 

a T Sov e tr si a eC f e °,^, StartS at the ,aser which ,s °" • "o»»t.tntop 3 km 

E ' velocity, VJ, can be adjusted In magnitude between 3 6 and 

10 km/sec by varying the energy ratio In thl two pulses 
F. The thruster efficiency will be ETH - 40* for any VJ In this ranoe 

th !‘ » Popellant can be found which will perfom a! *11 
«s? 1mplct;“of7ith?™. CalCU,ati ° n a,0l<Mn * environmental and 

G * I be a 5 i . re f t 'l 0 " thrust, which is normal to the vehicle base, can 

attit J H Sted J 1 ting the vehicle. It was assumed that the vehicle 
attitude would be continuously measured from the ground and controlled 

that Jf 9 th he J- e ! r PUlS6S ° ff center * A short simulation indicated 
that, if the vehicle was spinning at a few rps, it would be possible 

1 . to control * he thrust axis to within about 5 degrees. 

H. Vehicle design was assumed to allow an angle of incidence between the 

laser and the base up to 65 degrees without exposing the payload to 
damaging laser radiation (see Fig. 1). S payioaa to 

It was assumed that the beam director mirror would be 10 meters in 
dia. allowing the 10 micron beam to be focussed on the 1 meter di’a 
vehicle base out to a range of about 800 km. This impl es pe Romance 

SLElMr ° m the dlffraction ""It. More work il needed to 
specify tolerances on optical performance. 

The program used to calculate the trajectory is to be found in Appendix 1 

* r ™ UH *.*»" ‘HP numerical resultsare^erlp 

flble 2. The ascent to orbit is divided into four phases. 

1. Phase 1 starts with a liftoff close to the laser and with a vehicle 
weight (propellant plus payload) of 120 kg which was close to the 
largest load the laser could lift with VJ=3.6 km/sec. VJ was varied 
to make the best compromise between gravity and drag so as to minimize 

ihat m ^Mt° SS f e i U k U ?! titude 9 ain * Phase 1 was terminated (some- 
what arbitrarily) when the acceleration reached 1 g. At the end of 
phase 1 the mass was 57 kg. T 

2. In phase 2, continuing the vertical ascent to 130 km, the vertical 
acceleration was maintained at 1 g and VJ varied from 5 to 10 km/sec 
the mass ended up at 38 kg and the vertical velocity was 1.45 km/sec*. 

3. At the beginning of the extraatmospheric acceleration, phase 3, the 
vehicle was tilted so that the angle of incidence of the laser on the 

wa^at^h^K W ^ S • rad l an : The thrust, which is normal to the base, 

was at the beginning of phase 3, 33° up from the horizontal. The 

vehicle was maintained at this angle of Incidence to the laser as it 
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Figure 2. The sample trajectory. Beginning at a 3-km mountaintop, the launch is divided 
phases and reaches a 411-km orbit after 502 s. The vehicle coordinates and the angle (not scalable) to 
the laser are shown in the lower graph, and the vehicle acceleration is shown in the upper graph. 
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accelerated horizontally until after about 360 sec into the flight the 

b ? t r en the 1aser beam and the vertical reached .5 radians 
Then the thrust became horizontal and afterward had a downward compo- 
nent which was continued until the vertical velocity was cancelled. 

4 ' ?? en ^ vertical velocity became negative near the end of acceleration 
the vehicle was tilted in phase 4 to maintain the vertical velocity 

velocity L reached. ended ^ 5eC °" dS fnt ° th<i fll9ht " he " ° rbUa1 

TABLE 2 




SAMPLE 

10 MW LAUNCH 


TIME MASS 

HOR.DIST. 

HEIGHT 

VJ PHASE 1&2 

VERT. VELOC. 


PHASE 1 




0 120 

0 

3 

3.6 

.085 

20 107.09 

0 

4.7 

3.648 

.091 

40 96.4 

0 

6.6 

3.958 

.099 

60 87.45 

0 

8.7 

4.268 

.11 

80 79.77 

0 

11.1 

4.573 

.124 

100 73.13 

0 

13.8 

4.862 

.144 

120 67.25 

0 

17 

5.111 

.173 

140 61.87 

0 

20.9 

5.261 

.218 

160 56.59 

0 

26 

5.155 

.3 

173 

PHASE 2 




180 50.62 

0 

33.6 

4.973 

.468 

200 46.17 

0 

45 

6.575 

.664 

220 43.53 

0 

60.4 

8.055 

.861 

240 41.52 

0 

79.7 

8.761 

1.058 

260 39.74 

0 

102.9 

9.192 

1 . 254 

280 38.11 

0 

130.1 

9.588 

1 .451 

280 

PHASE 3 

HORIZONTAL VELOCITY 


300 35.86 

4.2 

159.8 

.412 

1 . 514 

320 33.61 

17.1 

190.7 

.863 

1 . 575 

340 31.36 

39.6 

222.7 

1.36 

1.619 

360 29.12 

72.7 

255.3 

1.924 

1.633 

380 26.89 

117.6 

287.8 

2.545 

1.605 

400 24.66 

175.6 

319.1 

3.231 

1 . 523 

430 22.46 

248 

348.2 

3.982 

1.377 

440 2028 

336.1 

373.5 

4.795 

1.157 

460 18.13 

441.1 

393.7 

5.672 

.855 

480 16.03 

564.3 

406.8 

6.608 

.465 

500 13.99 

706.8 

411.2 

7.605 

.02 

500 

PHASE 4 




502 13.79 

722.2 

411.2 

7.71 

.003 


END 

RESULTS 



INITIAL MASS, KG 120 


FINAL MASS 13.79 


RANGE= 831 

FINAL ZENITH ANGLE = 60 ACC. = 5 76 


ELEC. BILL/KG IN LEO = 

$ 10.1 

PROPELLANT = $ 15.39 
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Figure 3. 
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Estimates from Eq. (5) of laser 
from Avco for C0 2 lasers and 


launching capability to LEO (-400 km). Costs are based on 
from Itek for adaptive optics. 


Laser power, MW 



ECONOMICS OF THE 10 MW LAUNCHER 


What can we say now of the costs of transportation to LEO by this small 
scale laser propulsion? y 1 


The electricity used in the sample trajectory was about 505 kw hrs per kq 
of payload. Even though this is more than 50 times the ideal energy re- 

«ln me "f;- the C °?f ? f 55 ,s f'ectrlcity will „„t dominate transportation 
costs. This parallels the situation in chemical rockets where the fuel 

thfn T e pn 5° n ? t d ?Tn™ [ or the calculations below and for Table 2, 
the U.S. Govt rate of $„02/kw hr was used. 


rtl l Jo £ fain 7n te prope ' lant c ?? ta a realistic propellant Is 

n^inJ ? f ??• The u mass of Propellant, about 7.7 pounds per pound of 
ran hfl d f 1 s . sma ] 1 enough so that it would be expected that a propellant 

no? hL!T C J - heap enough so that propellant costs probably will 
not have an important impact on overall launch costs. For the calcula- 

remembe rid “t h» t l? ^ 2 ’ * 2/kg ° f P ro P e " a "* was assumed li mit be 
remembered that the choice of propellant will have a large impact on the 

oJ r a US f a % r pr f rd C1enC ^ and thUS - a dl>eCt impact on the 1 aunchi ng^capabi 1 i ty 
of a laser and on the economics of laser propulsion. a K y 

cost ;j f or laser propulsion are the capital and the operating 
H cos s . V he g: ound laser installation. Eq. 4 can be used to optimize the 
C ° S ^ between mirrors and lasers and to provided rough 
Mn^ramc f T J h ? capital . C0St of a la ser launching installation to launch 
ontir^tA K If f ° r e ? ampl ? we take the cost of mirrors made with adaptive 
T fill M h w P r °P°^ tlonal ^ the mirror area ($1 M/m- 2 was suggested' by 
Itek) and we take laser costs to be proportional to average (not peak) 

power, we get that the costs should be distributed equally between laser 
and mirror. If also we take a laser cost $25M + $ 5 per wattTestimated 

° f AVC ° f ° r C ° 2 lasers) ’ and correcting eq 4 by a 
factor 13.79/18 to agree with the result of the sample trajectory, we get 

Mo = 1 1 . 5*(C - 50)*C A . 5 ( 5 ) 

where C is the capital cost in millions of dollars. 

Eq. 5 is plotted in fig. 3. The Department of Energy uses a rule of thumb 

20* of thl ope £ atln 9 costs of large experimental installations of 

rnc+c* h c capital cost per year- If we add amortization of the capital 

oer^pIU 5 Thr?n E5 n - th ? ° f th ? gr0Und insta H a tion comes to 40% 

per year. The 10 MW installation would have a capital cost of $150 M and 
an operating cost of $60 M/yr. * 


^°?n5 he sam P 1a calculation, a 10 MW laser could launch 13.79 kg (30.4 lbs) 
in 502 secs. If the laser were used with a duty factor of 1 (62,821 launch- 
es/yr), it would then launch 866 tonnes/yr.Al lotting the $60M/yr costs 
to the payload launched gives for the 10 MW laser /y 


Launch cost/lb = $32/duty factor 
+ $12 (Electricity & Propellant) (6) 


The $1 000/lb, estimated (Ref. 6) for the mid 1990s ^®; ic *| l /o^’ laijnrh 
would be bettered if the duty factory were greater than .032 (2000 launch- 
es/yr). The break-even point in eq. 6 will change quite rapidly with laser 
power. Neglecting the favorable variation of the launch time and the elec- 
tricity and propellant costs with laser power, eq. 5 gives that a 20MW 
laser installation costing $250 M would launch 36 kg. The break-even 
would then°occur at a duty factor of .02 (1300 launches/yr). 

The primary uncertainty in these estimates comes from a lack of knowledge 
of what can be done to produce an efficient thruster without introducing 

too much flight hardware which has added so much *5® c °j£ ?I„ C ,521 Cal 
rnrkpts In the same trajectory it was assumed that the 40% thruster 
efficiency would be maintained down to a flux of 10 MW/sq cm. Propellants 
will need to be developed to achieve high thruster efficiency at low flux 
to make laser propulsion a serious contender for space transportation to 
LEO. In view of the fact that almost no effort has been devoted to this 
requirement it should be evident that a great opportunity exists to 
creatively design materials. 

It is a pleasure to acknowledge stimulating discussions of ^ is . 

with participants at the Livermore Workshop, especially Jordin Kare, Dennis 
Reilly and Rod Hyde. I am indebted to Freeman Dyson and Lowell *f°° d 
the important suggestion that primary emphasis be placed on finding the 
minimum system for an initial trial of laser propulsion to orbit. 
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Appendix 

A Sample Trajectory Using a 10-MW Laser 


ILASERPROP TRAJECTORY (USING THRUSTER EFF THEFF * 4) 

! USE 10 MW LASER 
! "LASPRP15" 12/31/86 

! INITIALIZATION 

OPEN #1:PRINTER 
OPEN #2:NAME "OUTFILE" 

ERASE #2 

LET M = 120E3 ! MASS 

LET MO* M 

LET CD * .4 ! ASSUMES DRAG LIKE A SPHERE 

LET B * 0 ! ACTIVATES PHASE 1 

LET AREA * 1E4 ! BASE AREA 

LET THEFF = .4 

LET Y = 3E5 ! ALTITUDE, MOUNTAIN 

LET YO = Y 

LET LPWR - lE14*(l-.r (1-EXP( (Y0-Y)/7E5) )/COS(TH) ) ! 10 MW, 10% VERT ATM LOSS 

LET VX - 0 ! HORIZONTAL VEL. 

LET RHO * (1.225E-3)*EXP(-Y/7E5) ! EXPONENTIAL ATMOSPHERE 

LET VY - SQR(2*M*983/(RHO*AREA*CD)) ! INITIAL VY TO MINIMIZE MDOT/VY 

LET VJ = 3.6E5 ! INITIAL JET VELOCITY, PHASE 1 

PRINT #1: , "TABLE 2, SAMPLE 10 MW LAUNCH" 

PRINT #1: 

PRINT #1: "TIME"; "MASS ", "HOR. DIST", "HEIGHT", "VJ PHASE 
1&2", "VERT. VELOC." 

PRINT #1: 

PRINT #1: , "PHASE 1" 

! TRAJECTORY 

FOR T = 0 TO 2000 

IF T/20 = INT (T/20) THEN 

PRINT #1: T; INT (M/10)/100;, 

PRINT #1: INT (X/1E4)/10, INT (Y/1E4)/10, 

IF B * 0 OR C * 0 THEN 

PRINT #1: INT (VJ/100)/1E3, 

ELSE 

PRINT #1: INT (VX/100)/1E3, 

END IF 

PRINT #1: INT (VY/1E2)/1000 
END IF 

IF T/5 - INT (T/5) THEN ! OUTPUT TO PLOTTER 

? PRINT #1: #2:INT (X/1E4)/10;",",INT (Y/1E4J/10 
! IF T>1 THEN PRINT #1: #2: T;",", VDOT 
END IF 

IF (VX A 2/(6.371 E8 + Y))>983*(1 + Y/6.371E8)A(-2) THEN ! ORBIT REACHED 

SOUND 500,1 

PRINT #1: T;INT (M/10)/100, 

PRINT #1: INT (X/1E4)/10,INT(Y/1E4)/10,!NT (VX/1E3)/100,INT (VY/1E2J/1000 
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EXIT FOR 
END IF 
» 

IF YC1.30E7 THEN ! VERTICAL ASCENT THRU THE ATMOS. 

LET A - 0 

LET RHO - (1.225E-3) *EXP (-Y/7E5) 

LET DRAG - CD* • 5*RHO*(VX A 2 + VY A 2)* AREA 
LET THRUST - (M'VDOT + M*983 + DRAG) 

LET VJ « 2* LPWR'THEFF/THRUST 
LET MDOT - THRUST/VJ 
! PHASE 1 USE LOW VJ FOR HIGH THRUST 
IF B - 0 THEN ! TO MINIMIZE MDOT/VY 

LET VDOT - SQR(2*M*983/ (RHO*AREA*CD)) - VY 
IF VDOT >983 THEN ! GO TO PHASE 2 

LET B - 1 

PRINT #1: T, "PHASE 2" 

EXIT IF 

END IF SET VDOT - 983 PHASE 2 

ELSE IF B - 1 THEN ! SET VDOT - 983 PHASE 2 

LET VDOT - 983 
END IF 

LET VY - VY + VDOT 

! OUT OF THE ATMOSPHERE PHASE 3 

ELSE IF VY>0 THEN ! TILT VEHICLE 1 RADIAN FOR HOR AND DOWN THRUST 

LET A - 1 
IP C * 0 THEN 

PRINT #1: T/'PHASE 3 HORIZONTAL VELOCITY" 

LET C - 1 
END IF 
LET VJ - 8E5 

LET MDOT - LPWR*2*THEFF/VJ A 2 

LET VDOT - MDOT*VJ/M 

LET VX - VX + VDOT*SIN(TH + 1) 

LET VY - VY + VDOT*COS (TH + 1J-983 + ( (VX)A2)/6.371E8 
ELSE IF VY<0 THEN ! PHASE 4 

IF D - 0 THEN 

PRINT #1: T, "PHASE 4" 

LET D » 1 
END IF 

LET A - 2 ! HOLD VY NEAR 0 

LET VJ - 8E5 

LET MDOT - LPWR*2*THEFF/VJA2 

LET VDOT - VJ*MDOT/M 

LET VYDOT - -983 + ( (VX)A2) /6.371E8 

LET VX - VX + VDOT*SQR(1-(VYDOT/VDOTOA2) 

LET VY - VY-.983 + ( (VX)A2)/6.371E8 + VYDOT 
END IF 

LET X - X + VX 
LET Y - Y + VY 
LET M - M - MDOT 

LET TH - ATN(X/Y) ! ZENITH ANGLE 

LET LPWR - 1E14*(1-.1*(1-EXP( (YO-YJ/7E5) )/COS(TH) ) ! 10MW, 10% VERT ATM 
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LOSS 
NEXT T 

LET D - IE-5*SQR(XA2 + YA2) ! RANGE IN KM 


PRINT #1: 

PRINT #1:" END RESULTS" 

PRINT #1: 

PRINT #1: "INITIAL MASS, KG"; INT (MO/1E3), 

PRINT #1: "FINAL MASS"; INT (M/10)/100 
LET AC * INT (.rVDOT)/100 

PRINT #1: "RANGE - ";INT(D),"FINAL ZENITH ANGLE 
PRINT #1: " ACC. - ";AC 


";INT(57.3*TH); 


LET EB - INT(100*T*5E4*.02/( (3600)*(M/1E3) ) )/100 
LET PB* INT(100*( (MO-M)/M*2)/100 ! $2/KG 

PRINT #1: "ELEC. BILL/KG IN LEO - $";EB 
PRINT #1: "PROPELLANT - $" PB 
END 


! 102/KWHR,LASER«20% EFF 
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